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Abstract. The dark-speckle method ( Labeyrie 1995 ) 
combines features of speckle interferometry and adaptive 
optics to provide images of faint circumstellar material. 
We present preliminary results of observations, and simu- 
lations concluding to the feasibility of exo-planet imag- 
ing from the ground. Laboratory simulations with an 
avalanche photodiode indicate the detectability of a stellar 
companion of relative intensity 10~^ at 5 Airy radii from 
the star. New, more general, expressions for the signal- 
to-noise ratio and integration time are given. Compar- 
isons with direct long-exposure imaging indicate that the 
method improves the detectability of circumstellar nebu- 
losity, faint companions and planets. 
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K*" ■ 1. Introduction 
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Extra-solar planets (Mayor &Queloz 1995, Marcy & But- 
ler 1996) can, in principle, be seen in ground-based images 
using the light cancellation in dark speckles to remove 
the halo of star light (Labeyrie 1995). Long-exposures 
with high-performance adaptive optics and a correction 
of the seeing-induced shadow pattern on the telescope's 
pupil are also proposed by Angel (1994). A less extreme 
atmospheric correction is needed with the dark-speckle 
method. Both methods are expected to reach the 10^^ 
relative intensity needed to detect a Jupiter-like planet 
near its parent star. Related space-based techniques are 



also considered (Malbet et al. 1995, Gezari et al 1997). In 
the longer term, resolved images of detected exo-plancts 
will in principle be obtainable even from the ground with 



long-baseline interferometric arrays (Labeyrie 1996). 
The " dark-speckle" method exploits the light cancellation 
effect occuring in random coherent fields according to the 
Bose-Einstein statistics. Although adaptive optics can re- 
construct the Airy peak, and possibly the first few rings, 
in the focal image of a bright star, the degree of "seeing" 



correction which it provides cannot be good enough to 
remove the fluctuating speckles in the surrou nding zone 



of the Airy pattern . Coronagraph i c devices (Lyot 1930 



Bonneau et al. 1975, Mauron 198C, Malbet 1996, Gay & 



Rabbia 1996, Roddier & Roddier 1997) can remove the 
steady and organized part of the straylight, i.e. the first 
few Airy rings if they emerge from the boiling speckled 
halo, and even if they are hurried but remain detectable 
with the dark speckle analysis. 

One is therefore left with the problem of extracting as 
much information as possible from this speckled halo. De- 
structive interference in the star light occasionally causes 
a dark speckle to appear here and there. When this hap- 
pens at the position of the planet's own faint Airy peak, 
the darkening cannot be as deep as elsewhere. The planet's 
Airy peak, also restored by the adaptive optics, indeed 
has a rather stable intensity, which adds to the star's lo- 
cal intensity. The intensity histogram is therefore locally 
distorted, and suitable algorithms can display the local 
distorsions in the form of a cleaned image. The exposures 
must be shorter than the turbulence life time. With a large 
telescope and hours of integration, a planet, typically 10^ 
times fainter than the parent star, is expected to become 
visible. 



2. Signal to noise ratio 

To calculate the method's sensitivity in a more general 
way than done by Labeyrie (1995), we first follow his 
derivation to the point where he calculates the signal-to- 
noise ratio. 

The different photon distribution from the star and the 
planet defines the signal-to-noise ratio SNR, and accord- 
ing to the central limit theorem: 



nP,(0)[l - Po(0)] = 5Ai?v/nn(0) 



(1) 
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where P<,{k^) and Po(fco) are the probabilities to detect fc* 
and fco photons originating respectively from the star and 
from the planet, per pixel in a short exposure, and where 
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n is the total number of short exposures. Replacing the 
values of P*(0) and i-b(O) in Eq. (1), we obtain: 



SNR = (1 



1 



(2) 



The calculation of signal to noise ratio given in Labeyrie 
1995 assumed /c* 3> 1, which is unnecessarily restrictive, 
and unrealistic in some of the cases of interest. As sug- 
gested by one of us (RR), a more general analysis can be 
made under the assumption that /cq <C 1 and /c* 3> k^. 
Equation (2) then becomes: 



SNR « fcr 



H- fc* 



if j is the number of pixels per speckles, thus: 



j + jh j + 



(3) 



(4) 



and GK^ 



Where Kq and K^, are the number of photons per speckle 
in a short exposure, respectively for the planet and the 

star. 

The variables used in Labeyrie 1995 were: 

Kq-G 

where T is the total integration time, t the short exposure 
time, R the star/planet brightness ratio, G the gain of the 
adaptive optics or the ratio between the Airy peak and 
the halo of speckles, iV* the total number of photons per 
second detected from the star. 
It provides a new expression for the SNR: 



SNR = Kq. 




J ^ G 



(5) 



The sampling j should be fine enough to exploit the dark- 
est parts of the dark speckles, for a given threshold of 
detection e, linking the performance of the adaptive op- 
tics (G) and the brightness ratio {R). The intensity across 
a dark speckle may be coarsely modelled as a cosine func- 
tion of the position p by: 



I{p) = cos 



(6) 



where Ih is the mean intensity and -v/J the speckle size. 
The intersection between I{p) and the line y — elh gives 
s, the size of the pixel over which the light is integrated. 
It limits the minimal measurable intensity elh (with < 

e <C 1). A detailed calculation gives s k, 

We can assume that e = G/ R. Indeed, if /q is the intensity 

e//j e//j Jq G 
of the Airy peak, e =^ = ^;^ = ;^ 



speckle size 




s : pixel size 

Fig. 1. Shape of a dark speckle according to Eq. (6). If It 
is the mean intensity, integration on a pixel of size s yields 
an intensity elh- 

We also assume in the following that a planet can be seen 
if its own intensity is higher than elh- 
Now we are able to calculate a value of j: 



J 



{X/DY 



[x/dY 



l.27^e{X/ DY 



R 

0.62- 



(7) 



Recent numerical simulations (Boccaletti 1998) have 
shown that R/G should not exceed 10^ to retain a rea- 
sonable value of the sampling parameter j. 
The recorded level is It, but with dark speckles this resid- 
ual level decreases to elh, where e depends mainly on the 
adaptive optics performance. 

With the value of j obtained, the SNR expression be- 
comes: 



SNR^^ 



tT 



R Y 0.62| + 
Solving for R, leads to a third-degree equation: 
GtTN? 



0.62i?3 + tN^R2 



[SNRY 



(8) 



(9) 



The Cardan method gives a single positive solution from 
which we derive a final expression for the integration time. 



fSNR\^ R , , 



(10) 



With the following values : R = 10^ D = Sm; G = 10^; 
SNR = 5; TO„ = 2.5; q = 0.2; AA = lOOnm and t = 20ms, 
we find 2.7 hours of integration, i.e. 50% more than the 
result given by Labeyrie. 

Angel's discussion (1994) of the long-exposure method 
leads him to the following expression: 



R 



G 



T 



SNR \ At 



(11) 
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where Atopt is the optimum short exposure time. Using 
again the same values, Eq.(ll) gives a limiting bright- 
ness ratio of 1.4.10^ and the typical brightness ratio of 
10^ would be reached in 140 hours. The long-exposure ap- 
proach would in principle be more sensitive if the adaptive 
optics and shadow pattern compensation could be made 
extremely good. It is however less sensitive with current 
levels of adaptive performance. 

In fact, the difference between the dark-speckle and direct- 
long-exposure methods is more subtle, and both work in 
different regimes. A critical value of the photon stellar flux 
{Nfc) can be easily calculated by equahzing Eq.(lO) and 
Eq.(ll) for the same value of the total integration time 
(T). It leads to a second degree equation in A^*, with a 
single positive root: 



G + VG^ + 4 X 0.62GR 
2t 



(12) 



Therefore, if the photon flux is above N^,c, the dark-speckle 
method is more efficient than the long exposure and, be- 
low this limit the direct imaging is better. N^c strongly 
depends on the adaptive optics. Let us take a numeri- 
cal exemple. If the goal is to reach a 10^ brightness ra- 
tio with a gain of 10^ and 20ms exposures, the photon 
flux must be higher than 1. 27. 10^ photons / s, which can 
be achieved with a large aperture and wide bandwidth 
(9.3.10^ph/s for = 0, D = 2.4m, AA = lOOnm and 
20% efficiency). However, for a space telescope with adap- 
tive optics, the speckle lifetime is under control with values 
of the order of Is for exemple. Thus, N^c is decreased to 
2.54.10^p/iotons/s, which is less restrictive. However, the 
telescope should not be so large as to provide a partially 
resolved image of the star, since it would fill-in the dark 
speckles. 



3. Dark speckle lifetime 

Because the number of photon-events detected per pixel 
in each exposure is critical, these exposures should be as 
long as possible without degrading too much the dark- 
ness of the dark speckles. The optimal value is obviously 
shorter than the usual speckle lifetime considered by Rod- 
dier, Gilli & Lund (1982). A tentative lower limit can be 
estimated by linearly scaling the speckle lifetime in pro- 
portion to the dark speckle size defined by Eq.(7). It leads 
to impossibly short exposures for large R values. However, 
since the adaptive optics decreases the wave disturbance, 
it increases markedly the speckle lifetime at positions close 
to the Airy peak, depending upon the factor n^, where n 
is the number of actuators across the pupil containing 



of them (Ryan 1996). The optimal exposure time is there- 
fore dependant upon the adaptive optics performance. 



Table 1. Brightness ratio and SNR obtained with a 
photon-counting avalanche photodiode, for two values of 
the sampling parameter j. The number of zero-photon 
events was counted on 250000 short exposures of 100/^s, 
totalling 25s of integration. The SNR was calculated from 
Eq.(l). 
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360000 
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9.3 






950000 


155959 


150017 


4.0 



4. Simulation and results 

To assess the dark-speckle method we did a laboratory 
simulation using a single-pixel photon-counting detector, 
in the form of an avalanche photodiode. The star and 
planet were simulated by two He-Ne lasers, with ad- 
justable attenuators. A Lyot-type coronagraph permitted 
to remove the star's Airy peak and rings, thus decreas- 
ing the local halo intensity 10 to 15 times. Artificial "see- 
ing" was generated with a moving scatterer, selected to 
provide a Strehl ratio approaching that typical of current 
adaptive optical systems. The equivalent peak/halo gain 
was G = 3.4.10^. The flux of the central star was 44.10*^ 
photons/s. Calculating an histogram of the detected pho- 
ton events, we determined the SNR by comparing the 
number of zero-photon events with the planet turned on 
and off. As listed in Table 1, the results strongly depend 
on the sampling parameter j. They are consistent with 
Eq.(7) which gives the maximum brightness ratio (Rmax) 
theoretically reachable. In these laboratory tests, the dark- 
speckle analysis outperforms the long exposure when the 
sampling exceeds 144 pixels/ speckle area. In this experi- 
ment the short exposure time (100/is) is about 100 times 
shorter than the speckle lifetime (10ms). Available photon 
counting camera do not yet allow quite as short exposure. 

We also used the CP40 photon-counting camera developed 
by Foy and Blazit (Blazit 1986). It has pixels of 50/Ltm 
and a rather low saturation level of 50000p/iotons/s. Each 
pixel has a lower dark noise than an avalanche photodi- 
ode, but a much slower response. We used the algorithm 
described in Labeyrie (1995), generating a "dark map" 
by counting, in each pixel, the number of 20ms exposures 
which contribute zero photon. As contributions from suc- 
cessive 20ms exposures are accumulated in the dark map, 
a planet's Airy peak is expected to emerge as a black dot 
among background noise. To obtain the cleaned image, re- 
sults are displayed in positive form using, for example, an 
inverse square law. 

The CP40 discriminates between events featuring zero 
photoelectron and those featuring one or more photo- 
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electron/pixel/exposure. Adding all exposures generates 
an image which brings out faint companions better than 
would a similarly long exposure on a CCD detector. 
The flux of the central star was 5.3.10^ photons/s, and 
the gain of the adaptive optics was about 1556. The op- 
tical system operated at A = 0.67 fim, F/D ~ 3200, and 
the mask diameter was 0.34", i.e. it covered the central 2 
rings of the Airy pattern. The planet was located near the 
fifth ring of the diffraction pattern. 

Figure 2 shows that a companion 966 times fainter is well 
detected in 116 seconds. The SNR measured on a speckle 
size region (37 x 37 pixels) is 799, while the dark-speckle 
model predict an SNR of 1008. This model does not take 
into account the halo shape which can explain the 20% 
discrepancy. These initial results, where the companion is 
brighter than the halo, are very modest with respect to 
the performance expected at a later stage, but have pro- 
vided useful insight for improving the instrument. 
Currently, as seen in Figs 2 and 3, the detection sensi- 
tivity is limited by the presence of spurious blobs in the 
cleaned image. These are caused by static aberrations and 
coronagraphic mask effects. These residual blobs may be 
substracted from data obtained on a reference star. 

Finally, dark-speckle data have been recorded at the 
152cm telescope of Haute-Provence, usin g, during a s ingle 
night, the BOA adaptive optics system ( Madec 1997 ) de- 
velopped by the Ofhce National d'Etudes et de Recherches 
Aerospatiales (ONERA). This system, optimized for vis- 



Fig. 3. Cleaned image, generated through dark-speckle 
analysis, of the multiple star HD144217(/3 Sco). A com- 
panion (/3 Sco B) appears near the masked star image, 
(mask diameter^ 0.5", F/D = 1200) 



companion 




Fig. 2. Laboratory simulation of stellar companion detec- 
tion with the dark speckle method. The artificial compan- 
ion, 966 times fainter than the star, is at the center of the 
dark map, emerging from the halo. A coronagraphic mask, 
which hides the star's central peak is visible in the lower 
right corner. The blob seen in the upper right part of the 
field is a "static speckle" caused by permanent aberrations 
and removable with a reference star. 



Fig. 4. Same image enhanced with unsharp mask and high 
pass filter to emphasize the contrast and median filter to 
smooth the image at the speckle scale. 

iblc light, reaches a Strehl ratio of about 0.4 at 0.6/im 
in long exposures and higher in short exposures. 30 min- 
utes of observation, with an interference filter centered 
at 0.67/iTO (AA — lOOA), evidenced the faint component 
of the spectroscopic binary HD144217 (a = 16/i05'26", 
S = -19°48'18", V = 2.62). On the detector, the fiux 
from the primary was only 11860 photons/s. The angu- 
lar separation is about 0.45" with an uncertainty due to 
the mask offset. Owing to the low elevation of the star, 
the adaptive optics gain was only 12. The SNR mea- 
sured on 16 X 16 pixels is 168 and allows to derive from 
Eq.(lO) a brightness ratio of 88, corresponding to a 4.8 
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magnitude difference. On the long exposure synthesized 
from the same data (1 photon-events analysis) the SNR 
is very similar, but a direct imaging should give a SNR 
of 42 according to Eq.(ll). However, a recent measure- 
ment of lunar-occultation (Evans 1983) gives Am = 3.3. 
In this case, dark-speckle analysis should provide a very 
good detection with a SNR of 756 instead of 168. Unfor- 
tunately, the Hipparcos mission failed to detect the com- 
panion, probably because the Hipparcos satellite is unable 
to achieve Am > 4 with small angular separation, which 
is consistent with our data. More observations are needed 
to verify the companion magnitude. 
A continuing observing program is initiated. 

5. Conclusion 

The simulations and tentative observations lead to the 
following remarks: 

1. One would like to sample as densely as possible to 
reach the bottom of the dark speckles, but there should be 
enough photons per pixel. The optimal sampling is there- 
fore critical and we guess that it should be about 500 pixels 
per speckle area. 

2. A fast photon-counting camera with a low dark 
noise, high saturation level and many pixels is needed. 

3. The observations required a narrow-band filter since 
the diffraction and speckle pattern are color-dependant. 
The speckles are themselves dispersed radially. To increase 
the bandwidth usable in speckle intcrferometry, Wynne 
designed a chromatic lens with magnification inversely 
proportional to wavelength (Wynne 1979). D. Kohler built 
a Wynne corrector and we found it efficiently applicable 
to the present situation, where the speckle's wavelength 
dependance is more nearly a linear scaling. The resulting 
smearing of the planet's peak is acceptable if the spectral 
band remains less than 100 nm. 

4. Different types of apodisation can be achieved, using 
a classical Lyot coronagraph, the interference coronagraph 
of Gay & Rabbia (1996), or the phase-mask coronagraph 
of Roddier & Roddier (1997). Both recent systems favor 
the detection of planets closer to the central star's Airy 
peak. Laboratory simulations with these varied devices are 
considered to compare their respective efficiencies. 

Our simulations verify the theoretical expressions given 
for the signal to noise ratio. The SNR measured from the 
photon- number variance (Eq.(l)), is consistent with the 
SNR expected from the model (Eq.(lO)). In these pre- 
liminary tests, we had to use an interference filter and low 
saturation level camera which provides a weak signal. We 
were consequently unable to reach enough sensitivity for 
detecting extrasolar-planets or even brown-dwarf compan- 
ions. 

The dark-speckle method is also applicable to space tele- 
scopes. Even without turbulence, optical defects create 
static speckles which can be made to fiuctuate with a few 
actuators, arranged in the form of an active optics system, 



or a fast random scatterer. We proposed a "dark-speckle 
camera" , the Faint Source Coronagraphic Camera for the 



Hubble Space Telescope (Gczari ct al 1997). The project 
is reconsidered for the New Generation Space Telescope. 
IR wavelengths are of interest for the detection of extra- 
solar planets, for two reasons: the planet's contrast is im- 
proved and, turbulence is easier to correct at these wave- 
lengths. The forthcoming developement of bidimensional 
sensors with low read noise should allow red and IR work. 
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